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We report a cylindrical-peristaltic shape transformation in axons exposed to a controlled osmotic 
perturbation. The peristaltic shape relaxes and the axon recovers its original geometry within 
minutes. We show that the shape instability depends critically on swelling rate and that volume 
and membrane area regulation are responsible for the shape relaxation. We propose that volume 
regulation occurs via leakage of ions driven by elastic pressure, and analyse the peristaltic shape 
dynamics taking into account the internal structure of the axon. The results obtained provide a 
framework for understanding peristaltic shape dynamics in nerve fibers occurring in vivo. 

PACS numbers: 47.20.Dr, 87.16.-b, 87.19.La 



Axons and dendrites, collectively known as neurites, 
are thin, tubular extensions produced by neuronal cells 
. Structurally, they consist of an outer lipid membrane 
sheath bound to a core made up of an elastic network of 
highly cross-linked biopolymers known as the cytoskele- 
ton. Neurites are known to lose their normal cylindrical 
geometry and become peristaltically modulated-a pro- 
cess commonly known as beading-under a wide range 
of situations. These include neuro-degenerative diseases 
like Alzheimer's 0, brain trauma Q, stretch injuries to 
nerves and in vitro as well as in vivo application of 
neurotoxins or drugs 0, among others. In stretch in- 
juries, tension in the axon is responsible for beading. In 
the other examples, disruption of cytoskeletal integrity 
appears to be a common feature. 

Peristaltic modes have been studied in tubular mem- 
branes and gels under tension [10,13. This "pearling 
instability" is driven by surface tension, as in the 
Rayleigh-Plateau instability of liquid columns "o*] . In liq- 
uid columns, peristaltic perturbations with wavelengths 
larger than the circumference grow as they reduce surface 
area at constant volume. Unlike liquid columns, cylin- 
drical vesicles/gels are unstable only beyond a critical 
tension, when the gain in surface energy due to a reduc- 
tion in area overcomes the elastic energy for deforming 
the membrane/gel. Similar arguments can explain peri- 
staltic modes observed in stretched nerves [T^ and cell 
protrusions exposed to toxins [ij ■ In all the above men- 
tioned biological examples the beaded state persists and 
no recovery has been observed. 

We report a novel, reversible, peristaltic instability ob- 
served in chick-embryo axons and in PC12 neurites 
after a sudden dilution of the external medium, as shown 
in Fig. The transition to the peristaltic state occurs 
at a threshold membrane tension, induced by osmotic 
swelling. However, unlike previous reports on beading, 
the osmotically induced instability relaxes and the neu- 
rite recovers its original shape within minutes. We show 
that the relaxation is due both to cell volume and mem- 
brane surface regulation. Remarkably, the shape stabil- 
ity depends critically on the swelling rate. We model 
the pearling/beading instability taking into account the 




FIG. 1: (a) Osmotically induced shape instability in chick- 
embryo neurons, (b) Image sequence showing the growth and 
relaxation of the instability in a PC12 neurite. 

gel-like dynamics of the cytoskeleton and the membrane- 
cytoskeleton coupling, and find wavelengths in agreement 
with experimental observations. Further, we propose 
that elastic stresses in the cytoskeleton provide the driv- 
ing force for water and ions to leak through channel pro- 
teins in the membrane |il| to effect volume regulation. 

Experiments.- The experiments are performed using 
a 10 X 5 X 1 mm'^ flow chamber with temperature main- 
tained within ±0.2 °C. The cells are grown on treated 
glass coverslips 12] and transfered to the flow-chamber. 
The flow velocity in the plane of the axon is about 1 
jims^^ and the shear rate is about 1 s^^, well below val- 
ues needed to cause any significant shear stress on the 
neurites. The flow is switched from normal medium, with 
an ion concentration (^o = 0.3 mol/1, to a diluted medium 
to induce a hyposmotic shock. Phase contrast images are 
recorded using a CCD camera at a magnification of 0.084 
/im/pixel. The neurite boundaries are traced with a typ- 
ical error of a couple of pixels using a home-developed 
image analysis software. The volume and area are com- 
puted assuming axial symmetry. Only neurites which 
are firmly adherent at the two ends and freely suspended 
along the length are used for measurements. 

After exposure to a hyposmotic solution the neurite 
immediately starts to swell with a uniformly increasing 
radius from its initial value Rq. Fig. |21 shows the typ- 
ical variation of the normalized volume V = V/Vq and 
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FIG. 2: Evolution of the normalized volume V and area ^ of a 
neurite with -Ro = 0.5 /xm at 33 °C. The ratio Vv/A (dotted 
hne) reflects any deviations from the cylindrical geometry. 
The inset shows the volume and shape responses when the 
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concentration is switched in the sequence (j>o — > 0.5<^o — > cf>o- 
The initial slow swelling is due to a non-steplike change in 
concentration, as revealed by using an absorbing dye. Note: 
the shape remains cylindrical for the 0.500 — > 4>o shock. 

area A — A/Aq. For sufficiently large osmotic shocks, a 
standing peristaltic mode sets in at a threshold radius Re- 
The plot of \/v I A shows the evolution of the peristaltic 

shape. An increase in \fv j A indicates a reduction in 
surface area compared to a cylinder with identical vol- 
ume. The mode amplitude increases with time, reaches 
a maximum and then decreases back to zero. The ampli- 
tude and the wavelength appear uniform along the entire 
length of the neurite (see Fig. . An analysis of the in- 
stability observed on more than a hundred PC12 neurites 
with lengths of ^ 50 — 500 /J,m reveals the following. 

(i) The instability occurs only above a critical dilution, 
which is about O.50o at 37 °C and about 0.7(/)o at 25 °C 
for an initial radius i?o — 0.7 /xm. The critical dilution 
decreases for smaller Kq or for lower temperature. 

(ii) The wavelength A increases with increasing i?o but is 
independent of the neurite length. For a given neurite, 
the wavelength remains approximately constant during 
both the growth and relaxation of the instability. 

(iii) The maximum volume Ym scales linearly with Vq, 
i.e., Kn/Vb is largely independent of the neurite selected. 

(iv) Typically, the instability sets in when the area in- 
creases by 3-8%. The maximum increase can reach 20%. 

(v) For osmotic shocks slightly below the critical value, 
no shape change occurs even when the area increases well 
beyond the 3-8% mentioned earlier. The growth rate of 
the area decreases with decreasing osmotic shock. There- 
fore, we conclude that the stability depends on the rate 
of area increase as well as its magnitude. In another set 
of experiments, when the dilution is performed in several 
steps of very weak osmotic shocks the neurites remain 
cylindrical even when the external solution is diluted to 
pure water. In this case the volume and area remain close 
to their initial values throughout the experiment. 

(vi) Water flow across the membrane is mainly through 
water specific channel proteins belonging to the aqua- 
porin family 13]. Blocking them using 500mM dimethyl- 
sulfoxide significantly reduced swelling and made the 



neurites more resistant to peristaltic shape change. 

Discussion.- Volume regulation. — We assume the 
neurite to be a cylindrical, isotropic gel with an adher- 
ing fiuid membrane on its surface. Based on fluorescence 
microscopy, we conclude that the membrane is firmly at- 
tached to the cytoskeleton at least during the onset of 
the instability. The osmotic pressure difference across 
the membrane is AH = kT((j)i — (jjg), where (j)e and 0; are 
the external and internal net ion concentrations. The 
volume flux of water through the membrane is 

= Lp (An - Ap) , (1) 

where Lp is the membrane permeability to water and Ap 
is the hydrodynamic pressure difference over the mem- 
brane. The gel response time is faster than the timescale 
of swelling, which is determined by Lp. The local nor- 
mal force balance at the membrane determines the hy- 
drodynamic pressure jump Ap = {K + 4^/3)Ai?/i?o + 
EAR/Rl = KAR/Ro, where R = Rq+AR is the neurite 
radius, K and /i are the compression and shear moduli, 
and E the stretching modulus of the membrane. The 
volume growth rate of the neurite is ^ = JwA. 

Volume relaxation occurs due to the opening of 
membrane-bound ion-channel proteins 0. This Regu- 
latory Volume Decrease is known to involve both "ac- 
tive" as well as "passive" processes, many of which are 
poorly understood and cell-type dependent 0|. In gen- 
eral, it is known that the ion conductivity G depends of 
the charge distribution across the membrane. However, it 
has been shown for PC12 cells that the membrane poten- 
tial vanishes during swelling induced by an hyposmotic 
shock and it does not recover during the entire relaxation 
phase 14, 15] . Hence, we assume a constant average G 
for simplicity and write the ion flux as 

Jion = Gfcrin^. (2) 

(Pi 

The time variation of the internal number N of ions is 
N = Jion^- Rescaling the conservation Eqs. (|1I2|) re- 
veals two characteristic timescales: = 2L^°n^> 

ter influx and rion = 2G{kTy ^'^^ ^'-"^ conduction, where 
H" « 700 kPa is the osmotic pressure of normal medium. 
For weak osmotic shocks the conservation Eqs. (|1I2I| can 
be linearized, resulting in a double exponential relaxation 
form for the volume, with characteristic timescales re- 
lated to and Tjon. At short times the volume change 
is given by AV = LpAUeAot, where AHg = AH(t = 0). 
Thus, the initial swelling is not influenced by elastic 
forces or ion flow. Lp measured from experiments in- 
creases from Lp = 0.8 - 1.6(±0.3) x IQ-^^ m/(Pa.s), in 
a 15-37 °C range. 

Numerical analysis of the volume response given by 
Eqs. H1I2() produces either a smaller swelling or a longer 
relaxation as compared to experiments. To account for 
these, we consider the following biological processes, (i) 
The open probability of mechano-sensitive ion channels 
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FIG. 3: Comparison between volume evolutions obtained 
experimentally and from the analysis, (a) For (j)o —> 0.7<^o 
shock. K = 0.025n°, Pa = 0.06n°, t„ = 22 s, and rion = 70 s 
for 25° C. K = O.OSn", Pa = 0.090°, r» = 18 s, and rjon = 60 
s for 36°C. (b) For 0.5<?!>o -» c?!>o, using K = 1.5n°, Pa = 0, 
Tiu = 34 s, and Tion = 24. The neurite is densely filled with 
cytoskeleton and organelles which have to be compressed close 
to their "dead volume" for V < Vo giving a higher effective 
value of K for (b). 

increases with membrane tension [J. Therefore, we ex- 
pect an increase in ion conductivity upon swelling. In- 
cluding a tension dependence on G, however, cannot ac- 
count for the fast relaxation and large swelling, (ii) Fil- 
amentous actin undergoes a transient depolymerization 
mainly due to an influx of Ca^^ followed by repolymeriza- 
tion and contractility, due to rnyosin activation by Ca^^ , 
during volume relaxation [llliil. Thus, this actin dy- 
namics may account for both increased swelling and fast 
relaxation. Indeed, including an active pressure Pa in 
Eq. ^ after the maximum in volume, gives quantitative 
agreement with experiments as in Fig. Ola). From the 
fits, we estimate G = 0.7 — 4 x 10^*/(Jm^s). Converting 
to electrical conductivity units gives G = 2 — 10 S/m^, 
in agreement with reported values 0|. We also esti- 
mate Pa ~ 0.05 — 0.111°. At sufficiently high tempera- 
tures, we often observe an "undershoot" of the volume 
(see Fig. which could be attributed to active con- 
traction of actin. This complex actin dynamics has not 
been reported for hyperosmotic shocks (sudden increase 
in external concentration). Indeed, the simpler analy- 
sis without active pressure fits well to our hyperosmotic 
shock data as shown in Fig. O^b). Shrinking responses 
have previously been analyzed for kidney cells in sus- 
pension assuming volume dependent conductivities and 
neglecting elastic pressure jlSf . 

Shape instability. — As the neurite swells, the mem- 
brane tension <t increases and for a critical swelling ra- 
dius Rc it is unstable to peristaltic shape modulations. 
In [m a scaling argument was developed to determine 
the fastest growing wavelength in cylindrical cells, as- 
suming that dissipation is due to viscous flow inside the 
cell. Here, we take into account the gel-like structure 
of the neurite cytoskeleton. The solvent flow inside the 
gel is given by — Vp -I- r^V^v — C(v — dtu) = 0, where 
V is the solvent velocity (V • v = 0), p the hydrody- 
namic pressure, u the gel elastic deformation, rj the sol- 
vent viscosity and ^ the gel-solvent friction. In addition 
the gel elastic force is balanced by hydrodynamic friction 
C(9tu - v) = ^V^u + {K + iM)V(V • u). The solvent 
flow outside the neurite is given by the Stokes equations. 
We consider the stability of a modulation in the neu- 
rite radius R = Rc + ee"*sin(gz), with amplitude e and 



growth rate (jj{q). We assume that the gel is in elas- 
tic equilibrium at the onset of the instability, i.e. radial 
compression modes have relaxed. 

The membrane imposes boundary conditions at the 
neurite surface. The solvent velocity is continuous, and 
the solvent flow through the membrane is governed lo- 
cally by Eq. (Q. The osmotic pressure difference is 
preserved within linear instability analysis, since a peri- 
staltic modulation conserves volume to linear order. Con- 
servation of membrane area implies dtUr+RcdzVz = j^. 
The total radial stress jump (hydrodynamic and elastic) 
at the neurite surface must balance the Laplace pressure 
due to membrane tension Acr^r = -PLapiace- Tangential 
force balance implies cr^ = £,{vz — dtUz), where is the 
gel elastic surface stress, and ^ the sliding friction coeffi- 
cient between the cortical actin network and membrane. 
Rescaling the equations reveals two timescales for the 
gel dynamics: the relaxation time of compression modes, 
Tc = CR^/K, and that of shear modes ts — Tj/fi. The 
compression modes relax slower than the shear modes, 
Tc oc [R/£)'^Ts, where £ is the mesh size of the gel. For a 
shear modulus fj, ^ 1 kPa and solvent viscosity 50 times 
water viscosity, ts ~ 10"'* s. 

We find that the membrane-cytoskeleton sliding fric- 
tion affects the instability to order £^/{(R). Assuming 
that the friction is due to anchored membrane proteins, 
and using common values for membrane viscosity, we es- 
timate that ^ <C CR, showing that ^ can be neglected. 

The neurite is unstable above a critical tension CTc = 
K+ifsfj. ^(9^)- Close to (Tc the fastest growing 

wavenumber is q — qRc — \J'^^^ with a corresponding 

growth rate w ^ ■^{ "'J'" )^, thus initially the dynam- 
ics is governed by the gel compression modes, which is 
qualitatively similar to the result in Above a second 
critical tension (Tc2 = GfJ-Rc, the instability grows with 
the much faster rate to ~ ^/\/tstc- The tension CTc2 
corresponds to the critical tension for pearling instabil- 
ity in an incompressible gel (ctc cTc2 when K ^ oo), 
which suggests that above ac2 the instability grows via 
a peristaltic shear deformation of the gel. For high ten- 
sion the fastest growing mode converges to the universal 
value q « 0.65, independently of the gel elasticity and 
gel-solvent friction [see Fig. E^a)] . 

Using K ^ 2Q kPa, estimated from volume relaxation 
analysis, and /i ^ 1 kPa 19] gives o'c2 ~ l-07crc and 
qc2 ~ 0.2. For i?o ~ 0.5 /im the critical tension for 
pearling is CTc ~ 3 x 10^** N/m. 

The above analysis sets an upper limit q < 0.65 on 
the fastest growing wavenumber which is in good agree- 
ment with observations [see Figs. 01a,b)]. Remarkably, 
we never observed wavenumbers below q < 0.2, even for 
dilutions close to the pearling threshold. This could be 
due to the steep increase in q close to tic- Also, the 
timescale for the growth of modes with small q, which oc- 
cur close to tJc, may be longer than the volume regulation 
time and hence these modes are not easily observable. 

Surface area regulation. — Figure |21 shows some re- 
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FIG. 4: (a) Theoretical plot for the fastest growing wavenum- 
ber with tension and histogram of observed wavenumbers. (b) 
Variation of A with Ro at 25 °C for (j>o — > 0.6700 ■ The scatter 
in the data is due to natural variations between cells. 

markable features: (i) the instability begins to relax be- 
fore the area and hence the apparent tension has reached 
its maximum value, (ii) the instability disappears when 
the measured area and volume are well above their val- 
ues at its onset and (iii) a 20% increase in area due to 
stretching alone is unexpected as most biological mem- 
branes are known to rupture beyond 3-5% of stretching. 
A reasonable explanation for these observations is that 
lipids are continuously added to the membrane thereby 
relaxing tension. Assuming that the membrane tension 
is the same at the onset and disappearance of the in- 
stability, we estimate a rate of increase of area of about 
4 X 10~^ /im^/s per unit length. This process would lead 
to a viscoelastic-like rate dependence of tension, giving 
rise to the critical rate of dilution mentioned earlier. The 
increase in area is most likely due to the fusion of tiny 
multi-lamellar vesicles to the outer membrane [23, l2l) . 
It is unlikely that this tension relaxation is due to flow 
of membrane from the extremities as (a) the onset and 
relaxation of the instability appears uniform along the 
length even in the longest neurites (~ 1 mm) and (b) we 
do not observe any flow of extraneous particles sticking 
to the outer membrane. 

Conclusion.- Our experiments on the dynamics of the 
pearling instability in neurites reveal the role of volume 
regulatory mechanisms and direct membrane tension reg- 



ulation in maintaining or recovering their normal cylin- 
drical shape. Recovery to cylindrical shape has not been 
observed in previous experiments on beading cells. A 
possible explanation is that osmotic swelling is a gentle 
perturbation to which the cell has time to react before the 
deformation results in irreversible damage. Indeed, we 
observe that the instability is irreversible for sufficiently 
strong osmotic shocks. Moreover, it is not just the mag- 
nitude of osmotic swelling which triggers the instability, 
but also the rate at which the cell swells. Thus, observa- 
tion of the instability is a direct way to probe the mem- 
brane surface regulation dynamics in neurites. We have 
also performed preliminary experiments on pulling neu- 
rites using a glass needle. There is a critical rate as well 
as magnitude of pulling beyond which peristaltic modes 
appear. For small amplitude modulations, the instability 
relaxes, whereas for strong modulations it persists. Phys- 
iologically, the volume and membrane surface regulation 
mechanisms should play an important role in stabilizing 
axon shape to sufficiently slow and gentle perturbations 
occurring in vivo. The theoretical analysis highlights the 
importance of cytoskeletal elasticity in the volume reg- 
ulation processes, and suggests that active contractions 
of the cytoskeleton may play an important role. The im- 
plications of our investigation to pathological conditions, 
especially to dendritic beading due to local inflammation 
[23, is an interesting topic for further research. 
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